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Abstract 

O. 

We systematically investigate the prospects of testing new physics with tau sensitive near de- 

X. 

tectors at neutrino oscillation facilities. For neutrino beams from pion decay, from the decay of 



radiative ions, as well as from the decays of muons in a storage ring at a neutrino factory, we dis- 
cuss which effective operators can lead to new physics effects. Furthermore, we discuss the present 
bounds on such operators set by other experimental data currently available. For operators with 



Q_i' two leptons and two quarks we present the first complete analysis including all relevant operators 

<d : 

simultaneously and performing a Markov Chain Monte Carlo fit to the data. We find that these 

effects can induce tau neutrino appearance probabilities as large as 0(1O -4 ), which are within 

>. 
\& reach of forthcoming experiments. We highlight to which kind of new physics a tau sensitive near 

in ' 

C"^ ■ detector would be most sensitive. 
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I. INTRODUCTION 

With the expected results of the Atlas and CMS experiments at the LHC, particle physics 
will enter a new era by directly exploring physics at TeV energies. Complementary to the 
direct tests, flavour experiments such as LHCb and SuperBelle will search for indirect signs 
of new physics. On similar time scales, neutrino physics is also entering a new era. Long 
baseline neutrino oscillation experiments are aiming at a first measurement of the remaining 
unknown leptonic mixing parameters 6*i3,<5mns, and the sign of Am| 1; and will determine 
the other parameters with unprecedented precision. 

For this purpose, neutrino oscillation facilities exploit beams from various sources: cur- 
rently from pion decay (conventional beams and super-Jbeams [l|, |2j) and in the future pos- 



sibly also from the decay of radiative ions (/9-beams 



or from the decays of muons in a 



storage ring at a neutrino factory [5|-l9(. In addition to measuring the leptonic parameters, 
neutrino oscillation facilities can also be excellent probes of new physics, especially in the 
lepton sector. 

To probe new physics, the near detectors of neutrino oscillation experiments are a powerful 
tool, in particular when t hey are capable of detecting tau leptons, as has for example been 



pointed out in Refs. [10|, |ll|. With the above mentioned neutrino beams, tau appearance 



(above a comparatively low background level) would signal new physics 



Already present beam experiments like MINOS [12[ could be suitable for performing 
such a new physics search with a tau sensitive near detector. This idea, known as the Main 

n 

Injector Non-Standard Interaction Search (MINSIS) [13] has recently been discussed, e.g., 



at Ref. 14|. The physics reach of such an experiment depends, from the theoretical side, on 
how well the relevant processes at neutrino production and detection are already constrained 
by other experimental data available. From the experimental side, it is mainly limited by 
the ability to discriminate a possible signal from the background. 

To explore the potential of testing new physics with near detectors at neutrino oscillation 
facilities, we systematically investigate the present bounds on higher dimensional operators 
that can produce tau appearance signals at such experiments. In Sec. [Ill we introduce the 
effective operator formalism used to encode the new physics effects that can lead to signals 
at near detectors and classify them in three types. In Sees. IHH HVl and [V] we discuss the 
present bounds on each of the types of operators and the extensions of the SM that can 



lead to the least constrained combination of operators. Finally in Sec. |VI]we summarize the 
present constraints on each of the processes and the probability level that a near detector 
would require to probe new physics in each sector. 

II. NEW PHYSICS AT NEAR DETECTORS 

In general, at a near detector experiment, non-standard neutrino interactions (NSI) can 
modify both the neutrino production mechanism as well as the neutrino detection. The NSI 
do not need to be flavour diagonal, thus leading to striking lepton flavour conversion signals 
at near detectors, where the standard neutrino oscillations have not developed yet. Indeed, 
v T could be produced in pion or muon decays via new flavour-changing effects, signaling 
physics beyond the Standard Model (SM) if they are measured at a near detector from 
such a beam. Conversely, neutrino detection via inverse /3 decay (neutrino-nucleon charged 
current scattering) can be affected and a v e or v^ can be detected in association with a r. 
It should be noted that, while the new physics we will discuss generally also induce NSI 
effects in the neutrino flavor propagation in matter, these effects do not have enough time to 
develop to be measurable in a near detector and we will therefore not discuss them further. 

In comparison to disappearance searches, which would be dominated by systematic errors 
on the neutrino flux normalization and cross sections as well as the response of the detector, 
the search for neutrino flavor appearance turns out to be an excellent probe of new physics. 
Moreover, the bounds on lepton flavour violation in the e-/i sector are generally stronger 
due to very well constrained processes such as muon to electron conversion in nuclei. In 
addition, the intrinsic SM backgrounds in the neutrino beams and the detector are much 
larger for the electron and muon flavours. Therefore, the prospects to search for new physics 
are best explored via neutrino flavour conversion effects to tau neutrinos. 

In order to parametrize the new physics that can give rise to a tau appearance signal we 
shall consider an effective theory approach where the Lagrangian has the form 

£ = £>SM + £ Weinberg + £d=6, (1) 

with Csm being the SM Lagrangian, ^Weinberg the d = 5 Weinberg operator giving rise to 



neutrino masses 



15( | , and Cd=6 a collection of effective gauge invariant d = 6 operators that 



can modify the neutrino production or detection in a way that leads to a tau signal at a near 



detector. There are several such operators 16|-|25|, depending on the production mechanism 



of the neutrino beam (tt decay for conventional neutrino beams and super-beam experiments, 
f3 decay for /3-beam and reactor experiments and \x decays for neutrino factories). We will 
classify the new physics operators in three different categories: 

• Four-fermion operators involving two leptons and two quarks. These operators, when 
involving a neutrino and a charged lepton, can modify the neutrino production through 
7r and (3 decay, as well as its detection via inverse /3 decay, thus leading to signals at 
near detectors. We will refer to these operators as 2L2Q. 

• Four-fermion operators involving only leptons. These operators can modify the neu- 
trino production through muon decays, thus leading to signals at near detectors in a 
neutrino factory facility. We will refer to these operators as 4L. 

• "Kinetic" operators involving a pair of Higgs doublets and a pair of leptons with a 
derivative. After the Higgs field develops its vacuum expectation value (vev) these 
operators contribute to the leptonic kinetic terms and, upon their canonical normal- 
ization, induce deviations from unitarity of the leptonic mixing matrix, modifying the 
W and Z couplings. Such a non-unitary matrix implies flavour conversion effects at 
zero distance since the flavour eigenstates are no longer orthogonal, thus leading to 
the desired signals at near detectors. We will refer to these operators as non-unitarity 
(NU) operators. 

Here we will study the constraints that the processes related by gauge invariance to the 
lepton flavour violating production and detection of neutrinos imply. Since neutrinos and 
charged leptons belong to the same SU(2) multiplets, lepton flavour violation (LFV) in 
the neutrino sector is usually accompanied by related processes in the charged lepton sector. 

l26l |27| that for the NU and AL operators, it is possible to 



However, it has been shown [17 



18 



cancel all contributions to processes with four charged fermions while lepton flavour violation 
is still present in the neutrino sector. For the 2L2Q operators, it is trivial to expand the 
operators in components and see that this cancellation is not possible for the four-charged- 
fermion (4CF) operators with quarks, since both up and down quarks carry electromagnetic 
charge and both contributions would need to cancel simultaneously. Therefore, in the study 
of these operators, care must be taken to include the bounds that are already present on 
semi-leptonic flavor violating processes. 



In the next sections we will use the bounds on lepton flavour violating processes to derive 
new bounds on the 2L2Q operators and review the existing bounds on the 4L and NU 
operators. We will also discuss the implications for near detector neutrino experiments, 
detailing the minimum oscillation probability level that would need to be probed in order to 
improve over the present constraints. We will also discuss briefly which possible new physics 
models could give rise to the least constrained operators described. 

III. NEW PHYSICS WITH QUARKS: 2L2Q OPERATORS 

A. Bounds on 2L2Q operators from r decays 

In order to derive bounds on the 2L2Q operators through flavour changing processes 
involving four charged fermions, we will use the following basis of gauge invariant opera- 
tors [I6j 

(ol Q )/ = [l"yl q ][q 1p q], (2) 

(Ol Q )J = [L"Yr a L a ][Q lp r a Q], (3) 

(Oed)J = [L"E a ][DQ], (4) 

(O eu )J = \L P E a ]ir 2 [QU] T , (5) 

where L a is the left-handed SU(2) lepton doublet with flavour a, Q is the first generation 
quark doublet while E a , D, and U are the right-handed charged lepton with flavour a, 
down quark, and up quark fields, respectively. The matrices r a refer to the generators of 
the SU(2) symmetry. Each of these operators O will be assumed to appear in Cd=& with the 
corresponding coefficient 2y/2GpC, i.e., 2v / 2Gf(Cec/)/* as the coefficient of (Oeu)^ anc ^ so 
on. Here, Gf is the Fermi coupling constant. Since the first two operators are self-conjugate, 
their coefficients are Hermitian, while the coefficients of the last two operators need not be 
(however, the conjugate term must also be added to the Lagrangian). 

In this paper, we are mainly concerned with the possible signals of new physics at a near 
tau detector at a neutrino oscillation experiment. The operator coefficients that contribute 
to these processes are 

( C Lq)/i T i ( C Lq)/' ( C E»)/i ( C Eu)J, ( C Ed)J^ and (^Bf/)/) 



for Vp beams, as well as the corresponding operators with the muons exchanged for electrons 
for v e beams, and we will study what bounds can be derived on them with present data. In 
order to compute these bounds, we consider the Particle Data Group (PDG) [28] bounds 
on several different lepton flavor violating decays of the r (see Tab. HJ 1 . Notice that there 
is some admixture of strange quarks in some of the mesons involved in the decays, e.g., 
in the i] meson. The bounds from such processes on the coefficients listed above should 
be correlated with the effective interactions with strange quarks. We therefore also include 
the possible NSI with strange quarks through the coefficients C\q,, C^q,, and Ces with the 
corresponding operators involving the second generation doublet Q' and the right-handed 
strange singlet S. In general, C LQ , and C L q, will only appear in the combination C L q, + C L q,, 
which is the coefficient of the operator combination that selects the strange component of 
the Q' doublet, so that the charmed quarks do not have to be taken into account. 

In order to compute the branching ratios, we need to compute matrix elements of the 
type 

M = (£U\C d=6 \r), (6) 

where II represents the meson (or mesons) involved in a given decay. For this purpose, we 
adopt the approach of Ref . 30[ , where the authors employed either the Partially Conserved 



Axial Current (PCAC) hypothesis 31 



-36 



, Vector Meson Dominance (VMD) model |37H3 



or chiral perturbation theory (%PT) 40H42| to the different Lorentz structures of quark 



currents present, in order to compute the hadron matrix elements. Resulting from this 
procedure are expressions for the partial r decay widths. These expressions consist of a 
prefactor, which is a function of measured quantities, such as the SM masses, and meson 
decay constants, as well as a combination of the different NSI coefficients relevant to each 
decay. The resulting decay widths are then compared to the SM decay width of the r 
in order to obtain theoretical predictions for the branching ratios, which are bounded by 
experiments. In Tab. [Ij we give the relevant operator combinations and prefactors appearing 
in each decay. We fully take into account the correlation among all the NSI coefficients. 

The general situation, where all coefficients are allowed to appear independent of each 
other, is rather involved. Thus, we will first give an example where we assume that Cex — 



1 Some of the bounds on the effective operators could be improved by // ± t t signal searches at hadron 
colliders in the future. See Ref. 



29( for a recent study. 



Process 



T -> tp 



T ^ £lu 



T -> 



r — J- &r 



£tj 



£K+K- 



£K°K° 



Prefactor 



1.7 



1.4 



0.84 



0.69 



0.20 



0.081 



0.014 



0.014 



Relevant combination of coefficients 



If 3 I 



If 1 I 



Ir 1 + r 3 I 



\ C LQ + ^ \ C ED ~Ceu]\ 



4 \ U ED ^EU 



t I 2 



■^+Clq - [Clq> +Clq>\ 



3m± f 1 , 

2-JT_ -jT [C £C/ + C £C ] - C BS 

Am s m T I 2 



3m 



'/ 



4m s m T 



J?2 



-jF 



C 1 



'E£> 






|C. 



BD 



-c 



££/ 



\^ED ^EU\ 



\c 



EU 



-ES\ 



■-EU 



'ES\ 



\Ced +Ces\* 



Y^ED ^^ES\ 



BR bound 



6.8 
6.3 



8.9 
1.1 



1.3 

7.3 



1.1 
8.0 



10~ 8 

io- 8 



io- 8 
io- 7 



IO" 7 

io- 8 



io- 7 
io- 8 



6.5-10" 8 
9.2-10" 8 



2.9 

1.2 



2.5 

1.4 



3.4 
2.2 



io- 7 
io- 7 



io- 7 
io- 7 



io- 6 

IO" 6 



TABLE I: Summary of the combination of2L2Q operator coefficients with their respective numerical 
prefactor in the r decay branching ratio. In the rightmost column, we quote the PDG bound on the 
process with the upper value corresponding to £ = fi and the lower to £ = e. Here, the parameters 

2 

are defined as uj t = ^ — m ^ rnj = 1.3 and — ^ = 1.4. The ratio of the decay constants are defined 



m T mu+m^ m s m T 

as &± = 3^5 and &' = *-jh^§S with F 8 = 0.154 GeV and F° = 0.025 GeV ]MJ, i.e., 
i^_ = 0.87, &+ = 1.4, and &' = 1.9. 



C EX with X E {U, D} and that all coefficients are real. In this scenario, the only parameters 
that appear are C\q and the combination Ced — Ceu- I n Fig- HI we show how different decays 
constrain this parameter space. Even if very simplified, this figure gives a flavour of the 
general situation, where the allowed regions are ellipsoids in the parameters on which they 
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FIG. 1: Allowed parameter region from the bounds on the shown r decays in the C\q-(Ced — Ceu) 
plane. It has been assumed that the coefficients are real and that Cex = C EX . 
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TABLE II: Bounds at the 90 % posterior probability on the operator coefficients from lepton flavor 
violating decays. The values have been derived through MCMC methods (see text for details). 



depend. The final bounds are then composed by combining the bounds from all considered 
decays. 

For practical purposes, we derive simultaneous bounds on all of the operator coefficients 
by adopting a Markov Chain Monte Carlo (MCMC) approach, using the MonteCUBES soft- 
ware [43] for the numerical simulations. This approach also allows the study of correlations 
among the parameters. The PDG bounds on the branching ratios are used as the width of 
our likelihood function, which we assume to be Gaussian and centred at vanishing branching 
ratios. We also assume flat priors for both the absolute values and arguments of the complex 
coefficients, as well as for the values of real coefficients. 

The derived bounds on the operator coefficients C after consideration of the full correlation 
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FIG. 2: The 68 ; 90 and 95 % posterior probability contours as an example of the correlations found 
among the parameters (Ced)// , (Ceu)// o- n d {Ced)h > see ^ ex ^ f or details. 

among them are given in Tab. [UJ The bounds on most of the parameters are essentially 
independent. However, there is a correlation between Ced and Ceu-, since in most branching 
ratios they appear in the Ced — Ceu combination, as can be seen in Tab. UJ This direction 
is thus more strongly bound as can be seen in the left panel of Fig. [2J A similar correlation 
is present for C ED and C EU . Conversely, since in the branching ratios Ced and Ceu tend to 
appear in the same combination as C ED and C EU , respectively, but in an incoherent sum, 
we find an anticorrelation between Ced or Ceu and C ED or C EU , since when one is small the 
present bounds allow the other to be larger. This situation is depicted for Ced and C ED in 
the right panel of Fig. |2j 

B. Implications for near detectors 

With the bounds on the effective operators derived in the former section, we will now 
discuss the expected number of signal events in a near detector experiment. Let us first 
focus on the case of a conventional v^ beam from n decay. The signal process is 

7T + -> (J, + l/ 

vN -»■ t~X, 



where N is a target nucleon in a detector and X represents all the final states of the neutrino- 
nucleon scattering. The u-r flavour change can enter at the neutrino beam source or in the 
detection process J44J |45[. The rate of this process is calculated with a coherent sum of the 



amplitudes 2 



AVtotai = (T-X\C C cWrN)(u T fi + \C d=Q \n + ) + (r' X\£ d=6 \^N) (z/ M /i + \C C c k + ) , (7) 



where the Lagrangian £cc is the charged current four- Fermi interaction induced by W boson 
exchange. The first term corresponds to the amplitude with the NSI at the source, and the 
second term to that at detection. Since the size of these NSI amplitudes is characterized by 
the coefficients C, which are constrained to (9(10 -4 ) (see Tab. [TTJ) , the ratio TZ between the 
r~ signal events and the lepton flavour conserving /i~ events induced by the SM interactions 
is naively expected to be of an order of the square of C, i.e., at most O(10~ 7 — 10~ 8 ). 
However, there is an enhancement mechanism of the NSI amplitude with a pseudo-scalar 
quark current in pion decays (see, e.g., Ref. 46(). This is due to spin conservation, since such 
an operator would involve fermions of the appropriate chirality and therefore not require a 
mass-suppressed chirality flip as in the SM. Indeed, the LFV pion decay rate would be given 
by 

2 



r(vr + -)► v T {T 



2{Cl Q ); + u, 



(Cl 



EDJV 



(.4 



EUIV 



r(7T+ ->• V^ 



(8) 



with the chiral enhancement factor uj^ defined as 3 



m„ 



m„ 



u ld 



rrta m u + m d 



21, 



(9) 



as derived from the PCAC [47]. From the MCMC global fit we obtain the bound 



{^-EDJIJ- \^E 



EUJV 



<4.2- 10" 



(10) 



for the relevant combination of coefficients. This translates into a constraint on the signal 
rate of 



K < 7.9 • 10" 



(11) 



2 There is a distance L between the neutrino beam source and the near detector in a real experiment, and the 
contribution from standard neutrino oscillation for v^ — > v T may not be negligible. The effect can be taken 
into account by coherently adding the oscillation amplitude (t~ X\Ccc\v T N)(v T \e~ lHL \vi l )(v l j jl u + \Ccc\' K+ ) 
into Eq. ([?])■ Here, H is the neutrino propagation Hamiltonian. 

3 The factor strongly depends on the value of the quark masses. Here we adopt the PDG average of 
m u + m d = 8.0 MeV taj. 



10 



The current direct bound on this process is given by the short baseline neutrino oscillation 



experiments NOMAD [48| and CHORUS [49(: 



K < 1.63 • 1(T 4 (NOMAD), (12) 

< 2.2 • 1(T 4 (CHORUS). (13) 

Recently, a new experimental proposal, MINSIS, for a tau near detector in the NUMI beam 



at Fermilab has been proposed [13j, [l4| . The expected sensitivity to 1Z is estimated to be 
O(l0~ 6 ) or even better. The sensitivity to the source and detection NSIs in the on-going and 
forthcoming conventional beam experiments with explicit implementation of near detectors 



was discussed in Ref . 50| . It shows that the sensitivity is limited with a near detector which 
does not have capability of tau-detection. 

The 2L2Q interactions would also affect the source and detection processes in a future 
/9-beam facility, where pure v e (or u e ) are produced through the decays of radioactive ions. 
With e-r flavour violating NSIs, the produced beam would be contaminated by v T . Thus, 
the process v e N — > t + X would signal these new physics at a near detector. The amplitude 
for the v T beam production can be written as 

AVsourcc = (^e + N'\C d=e \N) 

= 2V2G F (1) {(Cl Q ) e T {u(e)~f°P L v(v T )]+ [(4X + (4^] [u(e)P L v(v T )]} 
-2v / 2G F (a J )(C2 Q ) e r [M(e)YP LW (z/ r )], (14) 

in the non-relativistic limit for the nucleon part. Here, (1) and (<Tj) are the nuclear matrix 
elements for the Fermi and the Gamow- Teller transitions, respectively. Contrary to the case 
of conventional beams, there is no enhancement mechanism in the amplitude. Similarly, no 
enhancements of the NSI effects would be present in the detection process, which, as for 
the conventional neutrino beam, would occur through inverse j3 decay. It therefore follows 
that the 2L2Q search at a near detector of a beta beam facility would be a demanding task 
given the smallness of the signal. The sensitivity reach at the beta beam experiment with a 



particular setup of a near detector was studied in Ref. 51] . The flux in reactor experiments 



also originates from a beta decay and its sensitivity to the NSIs was discussed in Ref. [52] 



11 



C. Connection to models at high energy scale 



According to the discussion in Sec. IIIIBl we conclude that the most relevant 2L2Q 
operators for near detector searches are Oed and Oeu- in the following we will discuss 
which possible high-energy physics could give rise to these operators in their low-energy 
effective theories. 

The operator Oed can be decomposed into the following fundamental interactions 
through mediator fields 



17 



19|: 



(Oed\ 



[L E a }[DQi], with a colour singlet scalar $(2* 



■1/2J 



-l[L^YD c }[Q c ,^ p E a }, with a colour triplet (3) vector V 2 (2 v +5/6 ) 



rT^.^un 



— \[L •y p Q i ][D'y p E a \, with a colour triplet (3) vector UiiVi 



2/3 ) 



where the symbol Xy for mediator fields indicates their SU(2)l representation X, U(l)y 
hypercharge Y, and Lorentz nature Z G {s,f} with s for a scalar and v for a vector. The 
situation for Oeu presents a similar structure. 

The scalar mediator $ has the same charges as the standard Higgs doublet. Thus, in 
a general two Higgs doublet model (THDM), the additional Higgs doublets could play the 
role of $ and mediate the NSI and LFV interactions. Note, however, that the doublet may 
be an "inert" Higgs field which does not get a vacuum expectation value. The other two 
mediators are typically referred to as "leptoquarks" . 

In an explicit model to realise the 2L2Q operators, possible additional constraints 
would have to be taken into account and it is an interesting question whether the model- 
independent bounds derived in this study can be saturated in such a model. 

Similar types of models, and their phenomenological constraints, have been studied in 



53 



-|55|. Higgs me- 



the literature. For example, there are studies on LFV in the THDM 

56M60J| . embedding of the 



diated LFV in the minimal supersymmetric standard model [53|, 

doublet scalar mediator into a R-parity violating supersymmetric model [6jJ, and models 

with leptoquarks [62]. 
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IV. NEW PHYSICS WITH LEPTONS: 4L OPERATORS 
A. Bounds on AL 

As discussed in Sec. HI1 in the case of 4L interactions, a d — 6 operator exists such 
that neutrino production and propagation processes can be affected in a flavour violating 
manner, while the corresponding charged lepton flavour violation is not induced. As a result, 
processes such as the LFV r decays studied in Sec. IIHAI cannot be exploited to derive bounds 
on this other kind of NSI. Nevertheless, other constraints can be derived. 

In the case of 4L interactions, there is a d = 6 operator that avoids 4CF flavour violating 
processes [171 Il8|: 

(O ll )J = 2(J? a iT 2 L,)(L e iT 2 L c P). (15) 

As such, constraining this operator in a model-independent way can be challenging. In 
particular, its most relevant effect is to modify the \x decay. This process is experimentally 
measured with exquisite precision and used to extract the value of the Fermi constant Gp . 
Indeed, the measured value of Gp through \i decays would be modified in presence of the 
operators of Eq. ( TT5l) with coefficients {CLL)a as: 

G, = G F (|1 + (C LL )/| 2 + \{C LL )J\ 2 + \{C LL )/\ 2 + \{C LL ) T T \ 2 ) . (16) 

Hence, in order to derive bounds on the coefficients of these operators, the value of Gf 
should be compared with other process from which it can be extracted but are not similarly 
affected. Such an alternative is offered by the quark sector through f3 and K decays. These 
decays are used to measure the values of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements V U( i and V us and used to test the unitarity relation, experimentally [28| : 

\V ud \ 2 + |K S | 2 = 1-000 ± 0.001 (17) 

at 90 % CL. In order to make these measurements, G^ is used as a measurement of the Fermi 
constant. Thus, barring cancellations between the first order of {Cll)^ and the second order 
of the rest of the coefficients that violate lepton flavour, the latter can be bounded to be 
\(Cll)<x I < 0.032. For the signal ratio 1Z, this bound implies 

U < 1.0 • 10 -3 (G F measurement). (18) 
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B. Connection to models at high energy scale 



In the case of leptonic NSI, the operator of Eq. (fT5|) can be mediated by a charged scalar 
singlet lji 63], which is often contained in models for radiative neutrino mass genera- 



tion 



64 



671 ]. In this simplest realization, only three couplings of the scalar to the lepton 



doublets are different from zero and stronger constraints than the ones described in Sec. IIVAI 
apply. These constraints were studied in Refs. 26|, |63] and are summarized in Tab. IIHI along 
with the weaker model-independent ones. The operator can in principle also be realized by 
a combination of the other types of mediator fields, e.g., Z'(1q) and W(3q), in which the 
couplings between leptons and the mediation fields must follow a specific relation to cancel 
4CF J27]. 



V. NEW PHYSICS WITH KINETIC OPERATORS: NU 



Bounds on NU 



As for the kinetic operators, the effective d = 6 operator with Higgs doublets 

+ A0\T 



(Ouvv)J = (M'K^^A 



(19) 



leads to non-canonical kinetic terms for neutrinos after electroweak symmetry breaking 
and, after their canonical normalization, induces non-unitarity in the lepton mixing ma- 



trix 



20 



21 



681 ] . Through this particular operator non-unitarity is induced in such a way 



that only neutrino interactions are affected and the rest of the SM Lagrangian remains un- 
changed. For this reason, it is often called Minimal Unitarity Violation (MUV). Since the 
couplings to the SM gauge bosons themselves are modified, they mediate NSI that affect 
the neutrino production, detection and propagation in any experiment in the same way, al- 
though the gauge interactions for charged leptons, which are measured with a high precision 
and consistent with the SM, are not affected. Thus, the bounds from processes such as the 
ones studied Sec. IIII Al cannot be applied to NU operators, since they originate from the 
SU(2)l breaking with the vev of Higgs doublet. Bounds on the NSI of this form can be 
found in the literature ljj, |26|, I69H72] . and are also summarized in Tab. IIIH 



14 



Beam (channel) 


2L2Q 


4L 


NU 


n (/i — > r) 


7.9 • 1(T 5 


n/a 


4.4 • 10~ 6 


(e -»■ r) 


< 10~ 6 


n/a 


1.0- 10~ 5 


/" (/" -> t) 


< 1(T 6 


1.0 -KT 3 (3.2 -1(T 5 ) 


4.4 • 10~ 6 


/x (e -> r) 


< 10~ 6 


1.0 -lO" 3 (3.2 -10" 5 ) 


1.0- 10~ 5 



TABLE III: Bounds at the 90 % CL on the probability of tau appearance at a near detector in a 
neutrino beam from f3 decay, tv decay or [i decays for the three types of new physics. Two different 
values, the bound to the effective four Fermi interactions and that with the assumption of the singlet 
scalar decomposition (in parenthesis), are shown in the column of leptonic NSI. 

B. Connection to models at high energy scale 



The operator in Eg. (19) is a typical imprint of the mixing between neutrinos and heavy 



SM singlet fermions 



69 



(see, e.g., Refs. [73 



dimensional models 



70| , such as right-handed neutrinos in TeV scale seesaw mechanism 



741 ) or Kaluza-Klein modes of bulk right-handed neutrinos in extra 



. Other attempts to exploit the 



26 



221 



jWSB in NSI were also made with 



75 



76| . However, the construction 



a dimension eight operator with Higgs doublets 

of a concrete model with the dimension eight operators is either rather contrived or reduces 

to the construction at d = 6 discussed earlier. 



VI. SUMMARY AND DISCUSSION 

Tau appearance at near detectors of neutrino oscillation facilities provides a very promis- 
ing window to new physics. In this paper we have systematically investigated which types 
of new physics can give rise to such a signal and discussed the bounds on the correspond- 
ing effective operators from the experimental data currently available. These bounds are 
summarised in Tab. II III for neutrino beams from pion decay (conventional neutrino beams 
and super-beams), from the decay of radiative ions (/3-beams) as well as from the decays 
of muons in a storage ring at a neutrino factory and for the three different types of new 
physics specified in Sec. [TT1 The resulting bounds on the tau appearance probability have 
to be compared to the estimated background level which has an intrinsic value of 0(1CT 6 ) 
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with presently discussed r detection technology [77] . 

Regarding new physics effects leading to tau appearance in beams from 7r decay, only 
the operators with two quarks and two leptons (2L2Q) and kinetic operators leading to 
non-unitarity (NU) of the leptonic mixing matrix are relevant. While bounds on the NU 



operators have been discussed already in Refs. [ll|, |26|, I69H72J and have here only been re- 
viewed, for operators with two leptons and two quarks (2L2Q operators) we have presented 
the first complete analysis including all relevant operators simultaneously and performing a 
Markov Chain Monte Carlo fit to the data. Taking the chiral enhancement of certain opera- 
tors in 7r decay into account, we found that current bounds only restrict the tau appearance 
probability to be smaller than 7.9- 10~ 5 , which means that an experiment like MINSIS would 
be about two orders of magnitude more sensitive to new physics than existing experiments. 
We have also discussed how extensions of the SM with leptoquarks or additional Higgs dou- 
blets could give rise to such 2L2Q operators. On the other hand, present bounds on the 
NU operators restrict the tau appearance probability to be smaller than 4.4 • 10~ 6 , so that 
observing this signal above the background level becomes much challenging. 

For neutrino beams from (3 and \x decay, the chiral enhancement is not present and 
therefore the effects from 2L2Q are already too constrained from the current experimental 
data to lead to a signal above the background level. On the other hand, for tau appearance 
caused by NU effects, the bounds on the appearance probability are at the level of 10 -5 
and future neutrino oscillation facilities might be significantly more sensitive 
bounds on effects from NU operators are only slightly stronger for facilities with beams 
from 7r decays so some improvement can be achieved there as well. Finally, we would like 
to remark that, without including specific new physics generating these operators, there are 
operators with four leptons (4L) which could cause tau appearance probabilities at the level 
of 10~ 3 . However, when the operator is generated in an explicit model of new physics (i.e., 



8l|. The 



by a singly charged scalar field), then much stronger bounds apply |26|, [63] . 

In summary, we found that, although bounds from the current experimental data are 
already quite strong for some possible new physics effects, there are some types (c.f. Tab. HID 
of new physics operators to which tau appearance searches, already envisioned at present 
neutrino oscillation facilities, are very sensitive. A tau sensitive near detector might be able 
to probe new physics at sensitivity levels up to two orders of magnitude better than current 
bounds. 



16 



Acknowledgments 

It is a pleasure to acknowledge very interesting discussions, which motivated this work, 
with Adam Para and with the participants of the workshop organized by Belen Gavela in 
Madrid around the MINSIS idea. The work was also encouraged by discussions with Yoshi- 
taka Kuno and Shinya Kanemura. We also want to acknowledge discussions with Marc Sher 
and Deirdre Black on the computation of the hadron matrix elements employed here. This 
work was supported by the European Community through the European Commission Marie 
Curie Actions Framework Programme 7 Intra-European Fellowship: Neutrino Evolution 
[M.B.] and by the DFG cluster of excellence "Origin and Structure of the Universe". 



[1 

[2 
[3] 
[4] 

[5 

[6 

[7 
[8 

[9 

[10 

[11 

[12 
[13 



Y. Itow et al. (The T2K) (2001), hep-ex/0106019. 
D. S. Ayres et al. (NOvA) (2004), hep-ex/0503053. 
P. Zucchelli, Phys. Lett. B532, 166 (2002). 
Beta-beam task group 



website: http://beta-beam.web.cern.ch/beta-beam/task/index.asp 



S. Geer, Phys. Rev. D57, 6989 (1998), hep-ph/9712290. 

A. De Rujula, M. B. Gavela, and P. Hernandez, Nucl. Phys. B547, 21 (1999), hep-ph/9811390. 

A. Cervera et al, Nucl. Phys. B579, 17 (2000), hep-ph/0002108. 

A. Bandyopadhyay et al. (ISS Physics Working Group), Rept. Prog. Phys. 72, 106201 (2009), 

0710.4947. 

International design study of the neutrino factory 



website: http://www.hep.ph.ic.ac.uk/ids/ 



T. Ota and J. Sato, Phys. Lett. B545, 367 (2002), hep-ph/0202145. 

S. Antusch, C. Biggio, E. Fernandez-Martinez, M. B. Gavela, and J. Lopez-Pavon, JHEP 10, 

084 (2006), hep-ph/0607020. 

D. G Michael et al. (MINOS), Phys. Rev. Lett. 97, 191801 (2006), hep-ex/0607088. 

A. Para (2009), Main Injector Non-Standard Interaction Search MINSIS 

website: http://www-off-axis.fnal.gov/MINSIS/, 



[14] B. Gavela et al. (2009), MINSIS workshop in Madird 



17 



website: https://www.ft.uam.es/workshops/neutrino/default.html. 
[15] S. Weinberg, Phys. Rev. Lett. 43, 1566 (1979). 
[16] W. Buchmuller and D. Wyler, Nucl. Phys. B268, 621 (1986). 
[17] S. Bergmann and Y. Grossman, Phys. Rev. D59, 093005 (1999), hep-ph/9809524. 
[18] S. Bergmann, Y. Grossman, and D. M. Pierce, Phys. Rev. D61, 053005 (2000), hep- 

ph/9909390. 
[19] S. Davidson, S. Forte, P. Gambino, N. Rius, and A. Strumia, JHEP 02, 037 (2002), hep- 

ph/0112302. 
[20] A. Broncano, M. B. Gavela, and E. E. Jenkins, Phys. Lett. B552, 177 (2003), hep-ph/0210271. 
[21] A. Broncano, M. B. Gavela, and E. E. Jenkins, Nucl. Phys. B672, 163 (2003), hep-ph/0307058. 
[22] A. Ibarra, E. Masso, and J. Redondo, Nucl. Phys. B715, 523 (2005), hep-ph/0410386. 
[23] A. Abada, C. Biggio, F. Bonnet, M. B. Gavela, and T. Hambye, JHEP 12, 061 (2007), 

0707.4058. 
[24] C. Biggio, M. Blennow, and E. Fernandez-Martinez, JHEP 08, 090 (2009), 0907.0097. 
[25] C. Biggio, M. Blennow, and E. Fernandez-Martinez, JHEP 03, 139 (2009), 0902.0607. 
[26] S. Antusch, J. P. Baumann, and E. Fernandez-Martinez, Nucl. Phys. B810, 369 (2009), 

0807.1003. 
[27] M. B. Gavela, D. Hernandez, T. Ota, and W. Winter, Phys. Rev. D79, 013007 (2009), 

0809.3451. 
[28] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008), 



website: http://pdg.lbl.gov/ 



[29] T. Han, I. Lewis, and M. Sher, JHEP 03, 090 (2010), 1001.0022. 

[30] D. Black, T. Han, H.-J. He, and M. Sher, Phys. Rev. D66, 053002 (2002), hep-ph/0206056. 

[31] M. L. Goldberger and S. B. Treiman, Phys. Rev. 110, 1178 (1958). 

[32] S. L. Adler, Phys. Rev. 137, B1022 (1965). 

[33] W. I. Weisberger, Phys. Rev. 143, 1302 (1966). 

[34] D. J. Gross, S. B. Treiman, and F. Wilczek, Phys. Rev. D19, 2188 (1979). 

[35] A. I. Vainshtein, V. I. Zakharov, and M. A. Shifman, JETP Lett. 22, 55 (1975). 

[36] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B120, 316 (1977). 

[37] J. J. Sakurai, Annals Phys. 11, 1 (1960). 

[38] N. M. Kroll, T. D. Lee, and B. Zumino, Phys. Rev. 157, 1376 (1967). 

18 



[39] H. B. O'Connell, B. C. Pearce, A. W. Thomas, and A. G. Williams, Prog. Part. Nucl. Phys. 

39, 201 (1997), hep-ph/9501251. 
[40] S. Weinberg, Physica A96, 327 (1979). 
[41] J. Gasser and H. Leutwyler, Ann. Phys. 158, 142 (1984). 
[42] J. Gasser and H. Leutwyler, Nucl. Phys. B250, 465 (1985). 

[43] M. Blennow and E. Fernandez-Martinez, Comput. Phys. Commun. 181, 227 (2010), 0903.3985. 
[44] Y. Grossman, Phys. Lett. B359, 141 (1995), hep-ph/9507344. 
[45] M. C. Gonzalez-Garcia, Y. Grossman, A. Gusso, and Y. Nir, Phys. Rev. D64, 096006 (2001), 

hep-ph/0105159. 
[46] P. Herczeg, Phys. Rev. D52, 3949 (1995). 
[47] J. F. Donoghue and L. F. Li, Phys. Rev. D19, 945 (1979). 
[48] P. Astier et al. (NOMAD), Nucl. Phys. B611, 3 (2001), hep-ex/0106102. 
[49] E. Eskut et al. (CHORUS), Nucl. Phys. B793, 326 (2008), 0710.3361. 
[50] J. Kopp, M. Lindner, T. Ota, and J. Sato, Phys. Rev. D77, 013007 (2008), arXiv:0708.0152 

[hep-ph] . 
[51] S. K. Agarwalla, S. Rakshit, and A. Raychaudhuri, Phys. Lett. B647, 380 (2007), hep- 

ph/0609252. 
[52] T. Ohlsson and H. Zhang, Phys. Lett. B671, 99 (2009), 0809.4835. 
[53] P. Paradisi, JHEP 02, 050 (2006), hep-ph/0508054. 

[54] S. Kanemura, T. Ota, and K. Tsumura, Phys. Rev. D73, 016006 (2006), hep-ph/0505191. 
[55] S. Davidson and G. Grenier (2010), 1001.0434. 

[56] K. S. Babu and C. Kolda, Phys. Rev. Lett. 89, 241802 (2002), hep-ph/0206310. 
[57] A. Dedes, J. R. Ellis, and M. Raidal, Phys. Lett. B549, 159 (2002), hep-ph/0209207. 
[58] A. Brignole and A. Rossi, Nucl. Phys. B701, 3 (2004), hep-ph/0404211. 
[59] E. Arganda, M. J. Herrero, and J. Portoles, JHEP 06, 079 (2008), 0803.2039. 
[60] M. J. Herrero, J. Portoles, and A. M. Rodriguez-Sanchez, Phys. Rev. D80, 015023 (2009), 

0903.5151. 
[61] R. Barbier et al., Phys. Rept. 420, 1 (2005), hep-ph/0406039. 
[62] W. Buchmuller, R. Ruckl, and D. Wyler, Phys. Lett. B191, 442 (1987). 
[63] F. Cuypers and S. Davidson, Eur. Phys. J. C2, 503 (1998), hep-ph/9609487. 
[64] A. Zee, Phys. Lett. B93, 389 (1980). 

19 



[65] A. Zee, Phys. Lett. B161, 141 (1985). 

[66] K. S. Babu, Phys. Lett. B203, 132 (1988). 

[67] E. Ma, Phys. Rev. Lett. 81, 1171 (1998), hep-ph/9805219. 

[68] A. De Gouvea, G. F. Giudice, A. Strumia, and K. Tobe, Nucl. Phys. B623, 395 (2002), 

hep-ph/0107156. 
[69] P. Langacker and D. London, Phys. Rev. D38, 907 (1988). 
[70] P. Langacker and D. London, Phys. Rev. D38, 886 (1988). 

[71] E. Nardi, E. Roulet, and D. Tommasini, Phys. Lett. B327, 319 (1994), hep-ph/9402224. 
[72] D. Tommasini, G. Barenboim, J. Bernabeu, and C. Jarlskog, Nucl. Phys. B444, 451 (1995), 

hep-ph/9503228. 
[73] R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D34, 1642 (1986). 
[74] G. C. Branco, W. Grimus, and L. Lavoura, Nucl. Phys. B312, 492 (1989). 
[75] Z. Berezhiani and A. Rossi, Phys. Lett. B535, 207 (2002), hep-ph/0111137. 
[76] S. Davidson, C. Pena-Garay, N. Rius, and A. Santamaria, JHEP 03, 011 (2003), hep- 

ph/0302093. 
[77] D. Autiero et al., Eur. Phys. J. C33, 243 (2004), hep-ph/0305185. 
[78] E. Fernandez-Martinez, M. B. Gavela, J. Lopez-Pavon, and O. Yasuda, Phys. Lett. B649, 

427 (2007), hep-ph/0703098. 
[79] S. Goswami and T. Ota, Phys. Rev. D78, 033012 (2008), 0802.1434. 
[80] S. Antusch, M. Blennow, E. Fernandez-Martinez, and J. Lopez-Pavon, Phys. Rev. D80, 033002 

(2009), 0903.3986. 
[81] D. Meloni, T. Ohlsson, W. Winter, and H. Zhang, JHEP 04, 041 (2010), 0912.2735. 



20 



